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On the purpose of studying the effect of long-range Coulomb-interaction in strongly 
correlated electronic systems we bring in as its representative the nearest-neighbor repulsion 
{y) in addition to the on-site repulsion {u) and shall investigate the possibility of the super- 
conducting transition of carriers doped into the charge-density wave (CDW) state expected 
for V > u/4 in 2-dimensional square lattice. We shall see that strongly correlated hopping 
processes of doped carriers make the systems superconducting. The favored superconducting 
phase is of extended s-wave symmetry, and Tc ~100K is shown to easily be attained near 
the half-filling. 
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§1. Introduction 



In studying strong correlation among electrons culminating in high Tc cuprates, heavy 
fermion systems or organic conductors, one of the most intricate problems in condensed- 
matter physics currently, we call for the Hubbard model and its derivatives like t-J model. 
It is true that we have learnt much from these models but it is not at all certain whether the 
learning survives, even qualitatively, when the neglected long-range Coulomb interactions 
are restored. This is mainly because these models are, in spite of being simple, not tractable 
analytically, even in perturbation theory. To say the least of it we still do not know what 
the ground state of the two-dimensional Heisenberg antiferromagnet is, the most important 
starting point toward real understanding of high Tc copper oxides as commonly believed. 

In viewing the status it is worthwhile putting the problem other way around, namely, 
taking into account the long-rangeness of Coulomb interaction, to some extent, from the 
beginning. The new "Problemstellung" is physically more fit, especially for cuprates, because 
of their ionic character. For this purpose we have proposed ^-^ a model in which the nearest 
neighbor Coulomb repulsion, is brought in as an representative of long-range Coulomb 
repulsion, in addition to the on-site repulsion, u and the nearest neighbor hopping integral, 
t: 

^= X] (4aCj> + S>^i,a)+wX]ni,T^J,i+^ ni,anj,a' (1 - 1) 

<i,j>a i <i,j>a,a' 

where cJ^(ci,o-) is the creation (annihilation) operator of an electron of spin a at site i, rii^a 
is the corresponding number operator with nj the total electron number at the site and 
< i,j > stands for a nearest neighbor combination. We shall call this model as the t-{u, v) 
model. With the model we have discussed the possibility of Cooper-pair formation in 

We thought it better to avoid calling it as "extended Hubbard model" as sometimes done, since the same 
calling is often used for ones, e.g., including the hopping term to the next nearest neighbors, t' , the one 
which may be called as (t,t')—u model instead. In the same token the Hubbard model itself might be called 
as t—u model. 
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2-dimensional square lattice when doped and v > tt/4 ^ As shown there, an advantage 
of the model is that the ground state is known, in contrast to the conventionally assumed 
case of antiferromagnetic (AF) uniform charge distribution {v < u/4), to be the static, 
nonmagnetic charge-density wave (CDW) state in which lattice sites are separated into A- 
sublattice with two electrons in a spin singlet pair and B-sublattice with no electron (Fig.l). 
The observation is confirmed later numerically by the quantum Monte Carlo simulation ^) 
and by the exact diagonalization method Indeed at the absolute zero the ground state 
at half filling can be identified as an antiferromagnetic spin-density wave (SDW) state for 
w = u/v > A while for < 4, as the CDW state as predicted. Further, for w < 4, the CDW 
phase is shown to be rather rigid against doping even up to quarter filling Recently 
it is also shown ^) that the true ground state at half filling is the CDW state modified by the 
condensed charge fluctuation due to the hopping, just like the ground states of spin ordered 
states renormalized by spin waves. The single particle as well as collective excitaions are 
found to be gapped with rather wide openings (about 8v — u for the former and 7v — u for 
the latter). Thus, as long as \t\ ^ v or u , the thermal fluctuations do not disturb the 
ground state at temperatures even as high as v/ks ■ 

It is true that the system described by the model for w < 4 is non magnetic, charge-density 
wave ordered and exhibiting a small Fermi surface when doped, not necessarily in accord 
with experimental observations up to now, thus seems not relevant to high cuprates, as 
it stands. But, as mentioned above, it has definite advantages, namely the well defined 
ground state being rigid against doping and thermal fluctuation, so that it is worth further 
theoretical investigation. 

As such we would like to pursue the line set by ref.l) and study in this paper supercon- 
ductivity of the model when doped. Since the model is particle-hole symmetric we shall only 
consider the hole doping of Nh = 2N — holes where ATg is the total electron number and 
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2N is the number of whole sites. In §.2 we present the effective Hamiltonian for the doped 
holes given in ref.l) based on the second-order perturbation theory. Elimination of on-site 
processes is carried out systematically by introducing pseudo-spin operators in terms of on- 
site-avoiding hole operators (construction of the pseudo-spin operator is done in Appendix 
A). As a result of the elimination the Bloch energy of a hole and the hole-hole interaction 
are modified. Pairing of holes is favored, especially near the CDW-AF phase boundary, for 
the spin-singlet and s-wave like state as in the ordinary BCS mechanism. 

Due to the symmetry character of the gap function of the point group -D4), the 
corresponding gap equation is brought into a 3 x 3 matrix form in §.3 (its elements will 
explicitly be given in Appendix B). We shall solve the gap equation iteratively and present 
the results for the transition temperature, T^, as a function of hole concentration, rih = 
Nh/2N for various values of tu: It shows a maximum around rih ~ 0.05 and can reach 
a value ksTc ~ J/10 near the CDW-AF phase boundary {w = 4) with J = t^/4u, the 
exchange integral adopted by the ordinary Hubbard model, thus can easily be ~ lOOK. 
We also calculate the superconducting correlation function and the specific heat. It is found 
that the superconducting correlation length is about 2.8ao where ao is the lattice constant 
of the original square lattice. 

Finally in §.4 we give conclusions and discussions. §.2 Possibility of Cooper Pair Formatit 

Based on ref.l) we shall study the possibility of superconductivity of doped carriers (holes) 
in this section. 

As noted there, when doped, holes occupy the A-sublattice sites (corresponding to Cu"'", 
in case of Cu02 system) while B-sublattice (corresponding to Cu"'"^, ditto) remains intact 
unless doped heavily. Hopping integrals and potentials of holes to the second order in are 
summarized in terms of = u/v and g = t^/v in the following: 
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(2-1) 



(2-2) 



A) Hopping integrals: (Fig.2a) There are two free hopping integrals 

Ti^-2g/{6-w) , 
T2 = Ti/2 , 

and also two for correlated hopping ^ 

T^ = -g{l/{5-w) + l/{6-w)}, 
T2 = -9/{5 - w) . 

See Fig. 2b) for a typical correlated hoppig process. 

B) Hole-hole potentials: The potential energy Vm for holes of m-th neighbor apart is 
Vi = -8£f{l/60 + 1/(5 - -u;) - 3/(6 - w) + 2/(7 - w)} for m = 1 , 
V2 = Vi/2 for m = 2 , (2-3) 

V 1/^ = for m > 3 . 



To the above we must add, at least formally, hopping integrals and potentials including 
on-site repulsion, namely, T^°'^\ T2°^^ and Vq. All of them are of order of u and should 
properly be eliminated, as will be seen. 



Taking all the above into account and performing Fourier-transform over A-sublattice 
sites we get Hamiltonian for the doped hole system in the following form: 

n = no + ni , (2 - 4a) 

Ho = ' (2 -4b) 

k,a 

Hi = {2N)-^ J2 Vh{k,k')bl^bk+,A'+,,a'h',a' , (2 -4c) 

k,k' ,q,i7,a' 



Correlated hopping mechanism first presented in refl) is analogious to the superexchange mechanism of the 
attractive interaction between two localized spins which leads the system to the antifferomagnetic ordering. 
In both mechanism an intermediate state with higher energy plays an essential role and produces effective 
attraction between two carriers. 
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where N is the number of A-sublattice sites. Further 6^ ^ {bk,a) denotes the creation (anni- 
hilation) operator of a hole defined by 



RiEA 



^i,—a 1 

(2-5) 



where Ci,_o- (c| _^) is the creation (anihilation) operator of an electron which appeared in (1- 
1) and A represents the set of all the A-sublattice sites. In (2-4b) e{k) is the tight-binding 
Bloch-energy relative to the band bottom 



e (k) = — 2Ti (3 — cos kxa — cos kya — cos kxa cos kyO) 

_W,^ , , , , . (2-6) 
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- (3 — cos kxtt — cos kytt — cos k^a cos kya) 



where a = \/2ao ( oq is the lattice constant of the original square-lattice) denotes the lattice 
constant of the A-sublattice and W stands for the band width defined by 

W = 8\Ti\ = 16g/{6-w) . (2-7) 

Notice that the band bottom lies at 6Ti relative to the energy of a localized single hole^-*. 
Finally the hole-hole interaction, Vh, in (2-4c) can be given by the aforementioned hopping 
integrals and potentials as 

Vhik, k') = V(k, k') + 2T(fe, k') , (2 - 8a) 

where V{k, k') and T{k, k') are respectively given by 

V{k, k') = Vo + 2yi{cos(fc^ - kx)a + cos{ky - ky)a} 

+AV2 cos{k'y. — kx)acos{ky — ky)a , 

T(fe, k') = 2f[°''^ (cos kxa + cos kya + {k ^ k')) 

+4T2°^\cos kxtt cos kytt + (k k')) 

+4:Ti{cos{kx — k'^)a cos kya + cos k'^a cos{ky — k'y)a + {k k')} 

-\-AT2 (cos kxa cos k'ya-\- cos kya cos k'^a) . 
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(2 - 86) 



(2 - 8c) 



The tilded hopping integrals are defined by subtracting from the correlated ones the corre- 
sponding free hopping integrals, namely, % —Ti — Ti, since the latter are already taken into 
account in defining the Bloch energy. Thus, for example, Ti = T2 = —g{lj{^ — — 1/(6 — 
w)}. 

In order to discuss the possibility of superconductivity in the doped carrier system it is 
convenient to introduce pairing operators defined by 

ta{k) = b-k,a{-icr2cra)aa'bk,a'/V2 , a = 0, 1, 2, 3 , (2 -9a) 

where ctq is the 2x2 unit matrix and Ua for a = 1, 2, 3 are conventional Pauli-matrices. 
Needless to say to{k) corresponds to the singlet pairing while t{k) governs the triplet pairing. 
Further they satisfy 

to{-k) = to{k) and t{-k) = -t{k) . (2-96) 

In terms of the pairing operators the reduced Hamiltonian corresponding to Hi can be 
written as 

ni ^Hisp + -HiTP , (2 -10a) 

UisP = (2iV)-i ^+(^' 4(^0 *o(fe) , (2 - 10b) 

k,k' 

HiTP =(2iV)-i ^V_(fe,A;') t^(fe') > (2 -10c) 

k,k' 

where Hisp and Hitp are the reduced interactions respectively for singlet and triplet pair- 
ings and 

V±{k, k') = (V(A;, k') ± V(fe, -k'))/2 . (2 - lOd) 

We now concern ourselves with the elimination of "on-site" terms. In the case of triplet 
pairing we are automatically free from them owing to the exclusion principle so that we can 
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use HiTP as it stands with V_(fc, k') of (2-lOd) given by 

V-{k,k') = V^P\k,k') = 2Fi(sinA;^ a sin /c^a + sin kyU sin /c^a) 

+4V2 (sin A;a;a cos fcj^a sin k'^a cos A;ya+ 

cos kxtt sin /c^a cos fc^a sin kytt) (2 — 11) 

+8T1 (sin k^a sin A;^a cos A;^a + cos k'^a sin /s^a sin kytt 

+ sin /cjjtt sin k'^a cos /c^a + cos k^a sin /c^a sin kytt) . 
For singlet pairing, on the other hand, we have to restrict the Hilbert-space to the states 
of no doubly occupied sites. By so doing we can discard the terms explicitly pertainig to 
on-site quantities like Vq, T^"^K Further we have to eliminate the on-site processes inherent 
in the Fourier-transformed operators. To do so we replace, in a conventinal manner, the 
hole operators by 

(2 - 12a) 




RieA 



where h\ ^ and hi^^j are the on-site-avoiding operators defined respectively by 

(2 - 126) 

where fii^a = ^h^^ = Ci^-(jc\ _^ denotes the number operator of holes at site i of spin a . 
Then TCq of (2-4b) should be replaced by 

Ho = (2 -13a) 

k,a 

and Hisp by 

Hisp = N-^ Y.^+{k,k')il{k)io{k') , (2-136) 

k,k' 

with io{k) defined by the on-site- avoiding operators through (2-9a) and with V+(fc, A;') given 
as 

V+ (fe, k') = V(^) (fe, k') + V^'^^ (fc, k') , (2 - 13c) 



where 

V*-*-* {k, k') = 2Vi (cos k^a cos k'^a + cos kya cos k'ytt) 
+AV2 cos kxa cos k'^a cos kya cos fc^a 
+8T1 (cos kxa cos /c^a cos k'ya + cos /c^a cos /c^a cos A;^a+ 

(2 - lU) 

cos A;a;a cos A;j/a cos fc^a + cos /cajO cos kya cos fc^a) 
+8T2 (cos k^a cos A;^a + cos kya cos /c^a) 
-2(6Ti+£(fe)+£(fe')) , 

and 

V^^^ (fe, fc') = 4V2 sin fcaja sin A;^a sin kya sin A;^a . (2 — 13e) 

Note that the last term on the r.h.s. of (2-13d) comes from the fact that in ehminating the 
on-site potential term 

k k' Ri&A Rj£A 

Vq should be replaced by the Bloch energy of two localized holes, 2 x 6T1, and similarly in 
eliminating the on-site hopping term the term with T^°'^^ (k) but not with T'(o'^) (fc) should 
be taken out, namely 

2Ar-i^T(-)(fe)4(fe) 5^to(fc0, 

k k' 

thus, in effect -2{e{k) + 6T1) is left^). 

It is rather clear from the definition that the pairing interaction V*-^^ is of p-wave character 
while V(*) and V^'^) are of s- and d-wave character, respectively. In Figs. 4 a) c) we 
illustrate the pairing interactions in the first Brillouin zone for w = u/v = 3.9. As seen from 
them V^*) develops rather a wide and deep valley around the center of the Brillouin zone, but 
the attractive valley of V^*^^ is off the center and shallow ( IV^'^V^*^'''* I ~ 10""^ )> so that the 
s-wave like pairing seems favored. This is mainly because the correlated hopping processes, 
which give rise to large attraction in V^^\ are absent in V^''). As for the p— wave type pairing 
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V*^^) can be as large as V^*) in magnitude, thanks to the correlated hopping integral, but 
the former exhibits alternating valleys and hills around the center of the Brillouin zone, as 
can be seen from Figs. 4c) so that the transition temperature higher than that of s-wave is 
unlikely. Thus we shall study in the following the possibility of superconductivity for the s- 
wave pairing alone. It is worth noting here that the attractive valley of V*^*) becomes the 
deeper and wider the closer w is to 4, the CDW-AF boundary: This is mainly due to the 
correlated hopping terms in (2-11) which are seen from (2-2) most attractive at = 4 for 
tu < 4. 

Now the reduced Hamiltonian to be investigated is 

n(^) = 7io + iV-'5^V(^)(fe,fc')4(fc)^o(fe'). (2 -14a) 

k,k' 

For a consistent treatment of the Hamiltonian, the pseudo-spin formalism ^) is most suited. 
Its application is, however, not immediate since the hole operators were modified as to avoid 
on-site double occupation. As we shall show in Appendix A we can construct appropriate 
pseudo-spin operators from the modified operators, by making use of local mean field ap- 
proximation for the commutators among the modified pairing operators, b-k,ibk,^, 6^ |^-fe j. 
and ha{k) = h], J)k,a > ^is below: 



r^'\k) = ^^^^^^ {nT(fc)+^i(-fc)-l} , 

= (^^^feI,T&Ui , (2-15) 

The pseudo-spin operators satisfy usual SU (2) commutation relations; 

[f(+)(fc),f(-)(fc')] = 24,fc'T(^nfc) , (2- 16a) 

[f(^)(fc),f(±)(fe')] = ±5fc,fe'f(±)(fe) . (2 - 166) 
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In terms of the pseudo-spin, the reduced Hamiltonian (2-14a) including the chemical poten- 
tial term can be rewritten as: 



7^^*) = Y,{e{k)-ilh){2f^^\k) + l) + N-'J2'^SA^^^''>^^^^ (2-146) 

A* 5 A' 

where e{k) , jih and Vgj^^(fc,fe') respectively are the Bloch energy, the chemical potential 
and the s-wave like pairing-interaction modified by the introduction of the on-site avoiding 
pseudo-spin and are given by 

s{k) = (1 — ^)s{k) = W{3 — cos k^a — cos kya — cos k^a cos kyo) , (2-17a) 

with W ^{1- if)W , (2-17b) 



and 



Viff{k,k') = V^^\k,k'). (2-19) 



We can then readily write down the gap equation at a temperature T: 

AW(fe,T) = -AT-i J2 Viffik,k') A('\k,'T)Q{k') , (2-20a) 

fc'ec< 

with e{k) = ^ tanh(^)|^^^ , (2-20b) 
for the gap defined by 

A(*) [k, T) = -N-^ ^eff < ^^"^ (^) > ' (2 - 20c) 

fc'ec< 

where < • • • > denotes the thermal average 

<•••>= — ^ -7^, 2-21) 
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and f*^^) = ^(f^"''^ + ^) . We have taken the phase convention that < f*^^)(fc) >= with 
f*^^) = — f^~^) . Further E{k) stands for the Bogohubov's quasi-particle energy 



^ (2 - 22) 

with ^(fe) = s{k) — fXh ■ 

In (2-20a) and (2-20c) stands for the region in the first BriUuoin zone where the 

modified Bloch energy is less than an appropriate cut-off. We shall deal with the cut-off in 
the next section. 



3. Properties of the Superconducting Phase 

Based on the formalism given in the preceding section we shall solve the gap equation 
to determine the transition temperature and the energy gap, then study the superconducting 
correlation and the specific heat. We shall omit the superscript (s) since we restrict ourselves 
to the s-wave like pairing alone. 

The Transiton Temperature and the Energy Gap: Without any loss of generality we 
can factorize the gap function as 

A(fe,T) = d(T)c(fc,T) . (3 -la) 

We note that d{T) depends only on temperature and c(fc, T) depends both on temperature 
and wave vector . The wave-vector part, c(fc,T) , represents the symmetry of the gap and 
directly refiects the lattice symmetry and the anisotropy of the hole-hole interaction. Then 
the gap equation (2-20a) can be reduced to 

c(fe,T) = -l Veff{k,k')c{k',T)@ik') , (3-2) 
fe'ec< 
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where Veff{k, k') is given by (2-19) and(2-13d). The transition temperature is determined 
by the condition 

d{Tc) = . (3 - 3a) 

From the hole-hole interaction (2-13d) the wave- vector dependent part, c{k, T) , should take 
the following form 

c{k,T) = co{T) + ci{T) {cos k^a + cos kytt) + C2{T) cos k^a cos kytt . (3 — 4) 

Before going into the detailed study of the gap equation , we should note that, although the 
irreducible representations appearing in the right hand side of (3-4), 1, cos kxa + cos kya, 
and cos k^a cos kya , are independent and othogonal to each other, we cannot decouple (3-2) 
into each of these representations due to the presence of the thermal factor, Q{k') . We can 
only bring the gap equation into the 3x3 matrix form as 

M(T) c(T) = , (3 - 5a) 

where 

/co(T)\ 

c(T)= ci(T) , (3-6) 

V2{T)) 

(P{T) Q{T) R{T)\ 
M{T) = S{T) T{T) U{T) . (3-7) 
\V{T) W{T) Z{T) I 

All the matrix elements appeared in M{T) are listed in Appendix B. 

It should be noted here that there is an arbitrariness in the factorization, (3-la): 
Indeed an overall factor for c(fe, T) may be factored out and be absorbed into d{T) . Let 
us define the gap function at the center of the first Brilloiun zone (F point) as 

Ar(T) = d{T){cQ{T) + 2ci(T) + c^iT)) . (3 - 8a) 

t The form of the gap consists of irreduceble representations of r4 of the point group D4: Indeed 
(3-4) corresponds to the superposition of the s— (1) and extended s— wave solutions ( cos fca^a-l-cos fc^a , 
cos k^a cos kya). There is no symmetry-related distinction among these representations.^^ 
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Now one way for the factorization is to factor out (cq + 2ci + C2) from c(fe,T) and use 
Ar(T) in place of d{T) . Instead we factor out cq from c(fc, T) , for the sake of numerical 
simplicity, and redefine relevant quantities as 



' D{T) = co{T)d{T), 

C{k,T) = c{k,T)/co{T) = 1 + Ci{T){cos k^a + cos kya) + C2{T) cos k^a cos kya , 
with Ci = Cj/co , 

A{k,T) = D{T)C{k,T) = D{T){1 + Ci{T){cos ka;a + cos kya) + C2{T) cos ka;a cos kya} . 
Then Ar(T) is given by 



The gap equation, (3-5a), remains the same formally, but due to the above mentioned 
arbitrariness in the vector (c) , should turn out to be 



although M contains Cj in a nonlinear manner through E{k) . The coefficients, Cj , are 
self-consistently determined by 



Needless to say the transition temperature is determined from (3-3), or Ar(Tc) = , in 
other words, by 



Ar(T) = D{T){1 + 2Ci(T) + C2{T)) . 



(3 - 86) 



det M(T) = , 



(3 - 5b) 




(3 - 5c) 



det M(re)|Ar(T,)=0 = 0, 



(3 - 10) 



for a given set of hole concentration, , and w . 



To solve the gap equation we resort to the iterative procedure: To begin with we solve 
(3-10) and determine the transition temperture, Tc . Then at a temperature lower than Tc 
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we put the zeroth order solution for Ci{T) as 

Ci°^(T) = Ci(Te)|Ap(T.)=0, 

(3-11) 

(T) = C2(Te)UKT.)=0, 

where Ci{Tc) are determined by (3-5c) by setting T = Tc and Ar (Tc) = . We next write 
down the first iteration for A(fc,T) as 

A^^) (fe, T) = D^^^ (T){1 + Cj°^ (T)(cos k^a + cos kya) + cf ^ (T) cos /c^^a cos kya} , (3 - 12) 

and solve (3-5b) for D^^^T) . By using the solution we calculate the r.h.s. of (3-5c) to 
obtain the first iteration to Ci{T) . By repeating the procedure for higher iterations until the 
desired accuracy is attained we can reach the solution for the gap function at an arbitrary 
temperature below Tc . 

We shall now present the numerical results for the transition temperature as well as the 
gap function. As for the energy cut-off, ic , we take the energy contour corresponding to the 
Fermi surface at the hole concentration of Uhc = 0.25 : This value is chosen because above 
which the rigidity of the CDW ground-state for the t — {u, v) model is shown demolished, 
as stated before for Cu02 system^) . Then we have 

£c = 0.58W^ . (3 - 13) 

We shall call this as the natural cut-off. 

Firstly in Fig. 5a) wc give Tc as a function of nu for sevaral values oiw . As seen 
there it increases as ^^/rlh as n/^ increases from 0, then develops maximum at nu = 0.054 
for w = 4.0; at Uh = 0.048 for w = 3.9; at Uh = 0.043 and Uh = 0.192 for w = 3.8, 
respectively. It vanishes at n^^ ~ Uhc around the cut-off, as \/^n^f~^ri^ . Further, as w 
approaches to 4, the curve of grows up in accordance with the observation stated before 
that the pairing interaction becomes most attractive at the CDW-AF phase boundary. The 
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maximum value of the transition temperature reaches ~ 10~^, where the couphng constant 
g = t^/v corresponds to the conventional exchange energy, i.e., 4£^/u, of the Hubbard- 
model near = 4 and can be of order of ~ 10"^K. Hence our mechanism easily explains the 
transition temperature as high as ~ 100 K. We have also shown in Fig. 5b) the anisotropy 
parameters, Ci and C2, at Tg as functions of for the same set of values for w . Notice 
that 

Ci ~ C2 , (3 - 14) 

for all the values of w considered. Further they are rather large at low hole concentra- 
tions and at the CDW-AF boundary, decrease as rih increases and become small at high 
concentrations. 

As is pointed out just in the above, the transition temperature is dependent on the 
choice of the cut-off. To illustrate this we give Tg in Fig.6a) and Cj in Fig.6b) as functions 
of Hh at w = 3.9 for several choices of the cut-off. The cut-off dependence can clearly be 
seen in Fig. 6a). Further, as seen in Fig. 6b) the anisotropy is conspicuous at a low cut-off: 
This is because the main source of the pairing interaction is the correlated hopping terms 
in (2-13d), which are inherently anisotropic and give strong attraction around the center of 
the first Brilloiun zone, the region most important for a smaller cut-off. 

Next we provide the solution to the gap equation at temperatures below Tc at the 
natural cut-off. To be specific we have chosen w to be 3.9. To obtain an overall view of 
the gap in k space we plot on Fig. 7 A(/e, T) at the absolute zero and Uh = 0.05 . As seen 
from the figure the gap function is very much alike to the sign-reversed Bloch energy. This 
is owing to the approximate relation (3-14) which is seen holding also at T < Tg as shown 
in Fig. 8a). Thanks to this relation we can approximately express the gap function (3-9) in 
terms of the modified Bloch energy (2-17a) as 

A(fe, T) = D{T){1 + 3Ci(T) - 4Ci(T)(^)} . (3 - 15a) 

W 
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In other words the gap function has a constant value on a Bloch contour. We shall thus 
define from (3-15a) the energy gap, say Af{T) , of our model as the gap on the Fermi 
surface, e{k) = jlh , which is explicitly given by 

A^(T) = D{T){1 + 3Ci(T) - 4Ci(r)(a)} = Ar(T){l - )} • (3 - 156) 

In Fig.Sb) we have plotted Ai?(T) as a function of temperature for several values of hole 
concentration. The results for the gaps reveal general features of the mean field approxi- 



mation, namely, all of them vanish as 1 — (^)^ at Tg and develop wide plateau at low 
temperatures. Further in accord with the behavior of Tc as a function of Uh shown in 
Fig. 5a) the gap becomes the widest at nn = 0.05 at which Tc is the largest for w = 3.9 . 
Further we have presented in Fig. 8c) A_f(0) and Ar(0) , together with A(0) , the gap func- 
tion at the absolute zero averaged over the 2-dimensional Fermi sea, as functions of n/^ . Also 
shown in Fig.Sd) are the ratios 2Ar(0)//cBTc, 2Af{0) /ksTc and 2A{0) /ksTc as functions 
of hole concentration for w = 3.9. We can see that the gap on the Fermi contour gives the 
ratio close to the BCS value, 3.52. 

From (3-15b) we may find the condition for the gap on the Fermi surface to vanish, 
namely 

A, _ l + 3C,(r) _3 1 3 

W 4Ci(T) 4^4Ci(r) 4- ^ ^ 

From the above it is clear that the lower bound for ^ exceeds the natural cut-off set by 
(3-13) so that in our model the gap function never vanishes on the Fermi contour. For a 
hole concentration of our interest, i.e., < n^c the Fermi surface is much smaller than 
the contour on which the gap vanishes. 

Superconducting Correlation: As is already clear so far the gap function shows strong 
anisotropy refiecting the lattice symmetry and the anisotropic pairing interaction. It is then 
of great interest to know how this anisotropy affects the superconducting correlation. The 
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superconducting correlation function for the s-wave like pairing can be calculated at T = 



where r = {x, y) denotes the relative distance between paired holes in the A-sublattice coor- 
dinates. Recall that the coordinates are rotated by 90 degree with respect to the coordinates 
of the original lattice, say {X, Y), so that the x — y-line in the former corresponds to the 
y-axis in the latter and the a;-axis, to the X = y-line. On Fig.9a) we have plotted the 
correlation function along the X = y-line, respectively for rih = 0.03, 0.05, and 0.20 at 
w —3.9. It indicates the superconducting correlation falls rather rapidly as distance in- 
creases and reveals a clear correlation length, ^ , almost independent of hole concentration, 
which is estimated to be about ^ ~ 2.8 from the plot of maxima of xs (See Fig.Qb)). We 
have also found that the superconducting correlation in other directions is quite similar to 
the above, thus the superconducting correlation extends about ^ ~ 2.8ao in every direction. 

Electronic Specific Heat: Electronic specific heat in the superconducting phase is given 
by the standard expression; 



We have evaluated this quantity as a function of temperture and shown the result for = 3.9 
and rih = 0.05 in Fig 10. As a comparison we have also plotted there the heat capacity 



by 




(3 - 16) 




(3 - 17) 



in the normal phase. Cat- The jump ^ at T — Tc becomes 1.67 for w — 3.9 and 

rih = 0.05, a value little larger than the BCS value, 1.43. 
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4. Conclusions and Discussions 

Based on the 2-dimensional t — {u, v) model we have shown that the static , nonmag- 
netic CDW state, consisting of alternate doubly occupied sites (A-sublattice) and empty 
sites (B-sublattice), is eligible for superconductivity with Tc as high as lO^K when the 
system is, 

1) doped near the half-filling and 

2) close to the phase boundary to AF. 

Main source of pairing interaction between doped carriers, say holes which necessarily 
reside on A-sublattice (in case electrons, on B-sublattice), comes from the correlated hopping 
caused by consecutive hoppings of the background electrons via a nearby empty site (e.g.. 
Fig. 2b)). We have paid special attention to elimination of the on-site processes by making 
use of the pseudo-spin formalism. In accord with the elimination, the Hamiltonian of the 
doped system has to be modified. We have shown that the effective hole-hole interaction 
especially enhances the s-wave type attraction and results in the highest superconducting- 
transition temperature at the CDW-AF phase boundary, w — 4. 

Throughout this paper we have assumed that the background CDW configuration remains 
rigid against hole doping for rih less than Uhc = 0.25, the critical dose set by the numerical 
study for Cu02 system ^\ If more holes are doped, the CDW configuration will make a 
transition to a metallic phase. 

As noted in the introduction the CDW ground-state is shown modified by the con- 
densed charge fluctuation: There is in effect a fraction of charges on the originally empty 
sites, although the CDW state remains stable even at temperatures as high as v/ks ■ The 
finding is at the half-filling, thus it is of due interest to know how doping affects the ground 
state and the interaction among doped carriers. By doing so we can shed light on the 
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strongly correlated electronic system starting from the localized configuration, a picture 
complementary to the conventional ones. 

Lastly a few remarks should be made whether our model is relevant to real supercon- 
ductors. CDW configurations similar to ours are known to occur in substances such as 
polyacetylene(lD), NbSe2(2D) and BaBiOs (3D). In these systems the CDW instability is 
driven by the the Fermi-surfacc instability. For example BaBiOs system is a semiconduc- 
tor with the 3-dimensional CDW-configuration Ba2Bi^+Bi^+06. If Bi is substituted by Pb 
the system shows superconductivity at the Pb concentration of 65% or higher. The behav- 
ior of this compound can be explained in terms of the modified t-{u,v) model, although 
the Cooper-pair formation seems mediated by the electron-phonon interaction As for 
cuprate high-Tc superconductors it is believed that the ground state is in some AF phase, 
even doped, and no clue is yet found for any static CDW configuration. If, however, w 
is very close to 4, where our mechanism is most efficient for high Tc , the spin fluctuation 
might well emerge and play some role in the ground state as well as in the superconductive 
phase, since the boundary to AF phase is immediate. This possibility is to be investigated 
as a future problem together with the charge fluctuation mentioned above. 
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Appendix A 



We shall construct the pseudo-spin operator in terms of the on-site avoiding creation and 
annihilation operators of (2-12b). First we define r- operators by 

'^-"'"^ (A-1) 
Then we calculate the commutation relation of r -operators, i.e., 

Ri,Rj,Ri,Rm£A 

After some algebra we have for the commutator on the right hand side above 

- SimSji{l - Sij){l - ni,|)(l - %|) 

+ SijSiibl^bm,^ + SjmSjibl^bj^^ - Sjmbl^bl,l^,fi3,^ 

- ^iibl^hlbUm,^ + ^jm^iii^ - Sij)bl^hlb],ibj,^ ■ 



To the above we apply the local mean-field approximation for the number operators of holes, 
i.e.. 



Further, since the system is nonmagnetic, we can neglect terms which make spins fiip to- 
gether with terms ending up with 0(1 /N) to the main terms when Fourier transformed. 
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Thus we obtain 



where 



[r(+)(fc),r(-)(feO] =25,fe,(l-^)TW(fe) , (A - 2) 



TW(fe) = ^(nT(fe)+ni(-fc)-l) , with n^{k)^bljk,a . (A -3) 

In a similar manner we calculate commutators of r^^) with t^^^ to get 



Therefore if we define the pseudo-spin operators, f, by 

r(^Hfc) = ^(Y^(^T(fe) + ^i(-fc)-i)= (Y^r^^Hfe) , (A -5) 

and 

then we can show from (A-2,4) that the pseudo-spin operators satisfy the SU{2) commuta- 
tion relations, (2-16a,b). 
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Appendix B 



We shall provide the matrix elements appeared in the equation (3-7). By inserting the 
effective interaction (2-13d) and the gap function (3~4) into the gap equation (3-2) we can 
obtain them as follow: 

P{T) = 1 + 4(1 - —fTi{3A{T) - B{T) - C{T)) , 

Q(T) = 4(1 - ^)2Ti(3S(T) - D{T) - E{T)) , 

R{T) = 4(1 - !^)2Ti(3C(r) - E{T) - F{T)) , 

S{T) = -(1 - ^f{ViB{T) + 4Ti(2C(T) + B{T)) + 4:TiA{T)} , 

Zi 

T{T) = 1 - (1 - — f{ViD{T) + 4fi{2E{T) + D{T)) + ATiB{T)} , 

U{T) = -(1 - ^)^{ViE(T) + 4fi(2F(T) + E{T)) + 4TiC(T)} , 
z 

V{T) = -(1 - —f{2ViC{T) + 8TiS(T) + 4TiA{T)} , 

W(T) = -(1 - ^f{2ViE{T) + 8fiD{T) + 4TiS(T)} , 
z 

Z(T) = 1 - (1 - ^f{2ViF(T) + 8fiE(T) + 4TiC(T)} , 
z 



where 



A{T) = 1 , 



B{T) = cos kxtt + cos kytt 



C (T) = cos kxa cos kya , 



D{T) = (cos kxa + cos kya^ 



E(T) = cos kxa cos kya{cos kxa + cos kya) , 



F{T) = cos^ kx cos^ ky 
In the above denotes 



1 ^ tanhi/3E(/c) 
iV ^ 2E(k) 

k€C< ^ ' 
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Figure Captions 

Fig. 1: The static CDW configuration at the half-fiUing having alternate doubly occupied 
sites (A-sublattice) and empty sites(B-sublattice). Arrows stand for electrons of 
corresponding direction, and ao and a = V2ao respectively are the lattice costants 
of the original square lattice and A- or B-sublattice. 

Fig. 2: a) Classification of the hopping integrals. The hole to hop is marked by a circle, while 
the doubly encircled hole is a spectator; b) An example of processes leading to the 
correlated hopping integral, Ti . The dark arrows are for background electrons and 
the hatched ones for holes. By consecutive hoppings of a background electron (No. 3 
in the figure) via a neighboring empty site, a hole (No.l) in the initial configuration 
( I i > ) is effectively moved to the final one ( | / > ) through the intermediate one 
( I m > ). 

Fig. 3: a)The contour of the Bloch energy of doped carriers, e{k), in the first Brillouin zone 
for w = 3.9; b) The chemical potential, /ih, of doped carriers for w = 3.9. Energy is 
in unit of the band width, W = 16^/(6 — w) = 7.6g. 

Fig. 4: The pairing interactions in the first Brillouin zone along = —k'^ and ky = —k'y 
foTw = 3.9 respectively: a) of s— wave type, V^^^k,k') ; b) of p— wave type, 
V*-^-*(fc, fc') and c) of d— wave type, V^'^\k^k'). Energy is in unit of the band 
width, W = 16g/ {6 — w) = 7.6g. Notice difference in scale between c) and a) or b). 
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-^^g- ^- a)The transition temperature, T^, devided hy g — ^ and b) the gap-anisotropy pa- 
rameters, Ci and C2, as functions of hole concentration, n/j, respectively for w =4.0, 
3.9 and 3.8. The energy cut-off corresponds to the chemical potential at the natural 
cut-off for the hole concentration, n^c = 0.25. 

Fig. 6: a) The transition temperature in unit of g and b) the gap-anisotropy parameters, 
at w — 3.9 for various values of the energy cut-off set by the chemical potential 
corresponding to the cut-off in hole concentration, Uc- 

Fig. 7: The gap function, A{k,T), devided by g in the first Brillouin zone at T = for 
w = 3.9, Uh = 0.05 and the natural cut-off. 

Fig. 8: a) The gap-anisotropy parameters and b)the gap on the Fermi contour, Af{T) 
devided by , as functions of temperature for n/i=0.03 ,0.05 and 0.20. Also shown 

are, c) ^^^^^ Ar(o) A(o) ^\ 2Ar (0) 2Af(o) respectively as 

functions of hole concentration, rih • AH of them are calculated ai w = ^ — 3.9 and 
the natural cut-off. 

Fig. 9: a) The superconducting correlation function and b) the logarithm of its local maxima 
along the line X=Y for n/j=0.03, 0.05 and 0.20 at w =3.9 

Fig. 10: The specific heat in the superconducting ( Cs ) and normal ( Cjv ) phases normalized 
by Cs{Tc) for w = 3.9 and rih = 0.05 at the natural cut-off. 
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